
## Estimate 

data <- read_dta(file.path(dir, "DataCode/Data/Constructed/LISS/Clean_LISS_Wide.dta"))
data$choice <- data$choice == 1

# Hour dummies
data <- dummy_cols(data, select_columns = "wh", remove_selected_columns = TRUE)

#
data_men <- subset(data, woman == 0)
data_women <- subset(data, woman == 1)

model <- mlogit(choice ~ -1 + fl + tc + mw + ln_hwage + wh_1 + wh_2 + tc:mw + fl:mw,
                data = data_women,
                shape = "long", # = c("nomem_encr", "experiment", "choicealt"),
                choice = "choice",
                na.rm = TRUE, 
                alt.levels = 1)

print(model)